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Abstract. Regina is a software package for studying 3-manifold triangula- 
tions and normal surfaces. It includes a graphical user interface and Python 
bindings, and also supports angle structures, census enumeration, combina- 
torial recognition of triangulations, and high-level functions such as 3-sphere 
recognition, unknot recognition and connected sum decomposition. 

This paper brings 3-manifold topologists up-to-date with Regina as it ap- 
pears today, and documents for the first time in the literature some of the key 
algorithms, heuristics and implementations that are central to Regina's perfor- 
mance. These include the all-important simplification heuristics, key choices 
of data structures and algorithms to alleviate bottlenecks in normal surface 
enumeration, modern implementations of 3-sphere recognition and connected 
sum decomposition, and more. We also give some historical background for 
the project, including the key role played by Rubinstein in its genesis 15 years 
ago, and discuss current directions for future development. 



1. Introduction 

Algorithmic problems run to the heart of low-dimensional topology. Promi- 
nent amongst these are decision problems (e.g., recognising the miknot, or testing 
whether two triangulated manifolds are homeomorphic) ; decomposition problems 
(e.g., decomposing a triangulated manifold into a connected sum of prime mani- 
folds); and recognition problems (e.g., "naming" the manifold described by a given 
triangulation). 

For 3-manifolds in particular, such problems typically have highly complex and 
expensive solutions — running times are often exponential or super-exponential, and 
implementations are often major endeavours developed over many years (if they 
exist at all). This is in contrast to dimension two, where many such problems are 
easily solved in small polynomial time, and dimensions > 4, where such problems 
can become provably undecidable HU |3D] . 

Regina [5l [M] is a software package for low-manifold topologists and knot theo- 
rists. Its main focus is on triangulated 3-manifolds, though it is now branching into 
both two and four dimensions also (see Section [7]). It provides powerful algorithms 
and heuristics to assist with decision, decomposition and recognition problems; 
more broadly, it includes a range of facilities to study and manipulate 3-manifold 
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triangulations. It offers rich support for normal surface theory, a niajor algorithmic 
framework that recurs throughout 3-manifold topology. 

One of Regina's most important uses is for experimentation: it can analyse a 
single example too large to study by hand, or millions of examples "in bulk" . Other 
uses include computer proofs^ such as the recent resolution of Thurston's old ques- 
tion of whether the Weber-Seifert dodecahedral space is non-Haken ijS;, and calcu- 
lation of previously-unknown topological invariants, such as recent improvements 
to crosscap numbers in knot tables [TTl [TB] . 

The remainder of this introduction gives a brief overview of Regina and some 
history behind its early development, back to its genesis with Rubinstein 15 years 
ago. Following this, we select some of Regina's most important features and study 
the modern heuristics, algorithms and implementations behind them. Some of 
these details are crucial for Regina's fast performance, but have not appeared in 
the literature to date. Overall, this paper offers a fresh update on the significant 
progress and performance enhancements in the eight years since Regina's last major 
report in the literature 5 . 

The specific features that we focus on are: 

• the simplification of triangulations, which plays a critical role in many of 
Regina's high-level algorithms; 

• the enumeration of normal and almost normal surfaces, including the high- 
level O-efficiency, 3-sphere recognition and connected sum decomposition 
algorithms; 

• combinatorial recognition, whereby Regina attempts to recognise a trian- 
gulation by studying its combinatorial structure; 

• the enumeration of angle structures, taut angle structures and veering struc- 
tures on triangulations. 

With respect to these features, there are interesting new developments in the 
pipeline. Examples include "breadth-first simplification" of triangulations, the tree 
traversal algorithm for enumerating vertex normal surfaces, and new O-efficiency 
and 3-sphere recognition algorithms based on linear programming that exhibit ex- 
perimental polynomial-like behaviour. The corresponding code is already written 
(it is now being tested and prepared for integration into Regina), and we outline 
these new developments in the relevant sections of this paper. 

To finish. Section |6] gives concrete examples of how Regina can be used for 
experimentation, including an introduction to Regina's in-built Python scripting, so 
that interested readers can use these as launching points for their own experiments. 
Section [7] concludes with additional future directions for Regina, including new 
work with Ryan Budney on triangulated ^-manifolds, much of which is already 
coded and running in the development repository. 

1.1. Overview of Regina. Regina is now 13 years old, with over 190 000 lines of 
source code. It is released under the GNU General Public License, and contribu- 
tions from the research community are welcome. It adheres to the following broad 
development principles, in order of precedence: 

(1) Correctness: Having correct output is critical, particularly for applications 
such as computer proofs. Regina is extremely conservative in this regard: 
for example, it will use arbitrary precision integer arithmetic if it cannot 
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be proven unnecessary, and the API documentation makes thorough use of 
preconditions and postconditions. 

(2) Generality: Algorithms operate in the broadest possible scenarios (within 
reason), and do not require preconditions that cannot be easily tested. For 
instance, unknot recognition runs correctly for both bounded and ideal 
triangulations, and even when the input triangulation is not known to be 
a knot complement (whereupon it generalises to solid torus recognition). 

(3) Speed: Because many of its algorithms run in exponential or super-exponen- 
tial time, speed is imperative. Regina makes use of sophisticated algorithms 
and heuristics that, whilst adhering to the constraints of correctness and 
generality, make it practical for real topological problems. 

Regina is multi-platform, and offers a drag-and-drop installer for MacOS, an 
MSI-based installer for Windows, and ready-made packages for several GNU/Linux 
distributions. It is thoroughly documented, and stores its data files in a compressed 
XML format. Regina provides three levels of user interface: 

• a full graphical user interface, based on the Qt framework [48] : 

• a scripting interface based on Python, which can interact with the graphical 
interface or be used a standalone Python module; 

• a programmers' interface offering native access to Regina's mathematical 
core through a 0++ shared library. 

There are facilities to help new users learn their way around, including an illus- 
trated users' handbook, context-sensitive "what's this?" help, and sample data files 
that can be opened through the File — > Open Example menu. 

Regina's core strengths are in working with triangulations, normal surfaces 
and angle structures. It only offers basic support for hyperbolic geometries, for 
which the software packages SnapPea [59] and its successor SnapPy [23] are more 
suitable. In addition to its range of low-level operations and 3-manifold invari- 
ants, Regina includes implementations of high-level decision and decomposition 
algorithms, including the only known full implementations of 3-sphere recogni- 
tion and connected sum decomposition. For a comprehensive list of features, see 
http : //regina. sourcef orge .net/docs/f eatureset .html 



1.2. History. In the 1990s, Jeff Weeks' software package SnapPea had opened up 
enormous possibilities for computer experimentation with hyperbolic 3-manifolds. 
In contrast, outside the hyperbolic world many fundamental 3-manifold algorithms 
remained purely theoretical, including high-profile algorithms such as 3-sphere 
recognition, connected sum decomposition, Haken's unknot recognition algorithm, 
and testing for incompressible surfaces. 

More broadly, normal surface theory — a ubiquitous technology in algorithmic 
3-manifold topology — had no concrete software implementation (or at least none 
that had been publicised). Many of the algorithms that used normal surfaces were 
so complex that it was believed that any attempt at implementation would be both 
painstaking to develop and impractically slow to run. 

Nevertheless, Rubinstein sought to move these algorithms from the theoretical 
realm to the practical, and in 1997 he assembled a team for a new software project 
that would make computational normal surface theory a reality. This initial team 
consisted of David Letscher, Richard Rannard, and myself. Much planning was 
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done and a little preliminary code was written, but then the team members became 
occupied with different projects. 

Letscher later resumed the project on his own, and at a 1999 workshop on com- 
putational 3-manifold topology at Oklahoma State University he presented Normal^ 
a Java-based application that could modify and simplify triangulations, and enu- 
merate vertex normal surfaces. This was an exciting development, but the software 
was proof-of-concept only, and later that month Letscher and I sat down to begin 
again from scratch. The successor would be a fast and robust software package 
with a G++ engine, would be designed for extensibility and portability, and would 
support community-contributed "add-ons" . 

Letscher stayed with the project for the first year, and then moved on to pursue 
different endeavours. I remained with the project (a PhD student at the time, with 
more time on my hands), and the project — now known as Regina — saw its first 
public release in December 2000. Development at this stage was taking place at 
Oklahoma State University, where Bus Jaco was both a strong supporter and a 
significant influence on the project in its formative years. 

Around that time, a significant phase shift was taking place in normal surface 
theory: Jaco and Rubinstein were developing their theory of 0-efficient triangula- 
tions. Their landmark 2003 paper |33| showed that key problems such as 3-sphere 
recognition and connected sum decomposition could be solved without the expen- 
sive and highly problematic operation of cutting along a normal surface, and that 
instead this could be replaced by a crushing operation that always reduces the 
number of tctrahedra. For the first time, the key algorithms of 3-sphere recognition 
and connected sum decomposition became both feasible and practical to imple- 
ment. Shortly after, in March 2004, Regina became the first software package to 
implement robust, conclusive solutions to these fundamental problems. 

Regina has come a long way since development began in 1999. It continues to 
grow in its mathematics, performance, usability, and technical infrastructure, and it 
enjoys contributions from a wide range of people (as noted in the following section). 
For a list of highlights over the years, the reader is invited to peruse the "abbreviated 
changelog" at http : / / regina . sourcef orge . net /docs /hi story . html, 

1.3. Acknowledgements. Regina is a mature project and has benefited from the 
expertise of many people and the resources of many organisations. 

Ryan Budney and William Pettersson have both made significant and ongoing 
contributions, and are both now primary developers for the project. A host of other 
people have contributed either code or expertise, including Bernard Blackham, 
Marc Culler, Dominique Devriese, Nathan Dunfield, Matthias Goerner, William 
Jaco, David Letscher, Craig Macintyre, Melih Ozlen, Hyam Rubinstein, Jonathan 
Shewchuk, and Stephan Tillmann. In addition, the open-source software libraries 
Normaliz [3] (by Winfried Bruns, Bogdan Ichim and Christof Soeger) and SnapPea 
[59] (by Jeff Weeks) have been grafted directly into Regina's calculation engine. 

Organisations that have contributed to the development of Regina include the 
American Institute of Mathematics, the Australian Research Council (Discovery 
Projects DP0208490, DP1094516 and DP110101104), the Institute for the Physics 
and Mathematics of the Universe (University of Tokyo), Oklahoma State University, 
the Queensland Cyber Infrastructure Foundation, RMIT University, the University 
of Melbourne, The University of Queensland, the University of Victoria (Canada), 
and the Victorian Partnership for Advanced Computing (Australia). 
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2. Simplifying triangulations 

The most basic object in Regina is a 3-manifold triangulation. Regina does 
not restrict itself to simplicial complexes; instead it uses generalised triangulations, 
a broader notion that can represent a rich array of 3-manifolds using very few 
tetrahedra, as described in Section \2A\ below. 

One of Regina's most important functions is simplification: retriangulating a 3- 
manifold to use very few tetrahedra (ideally, the fewest possible). Regina includes 
a rich array of local simplification moves, as outlined in Section 12.21 together it 
combines these moves into a strong simplification algorithm, which we describe in 
full in Section [2731 Moreover, this algorithm is fast: in Section [2741 we outline some 
important implementation details, and show that it runs in small polynomial time. 

This simplification algorithm gives no guarantee of achieving a minimal trian- 
gulation (where the number of tetrahedra is the smallest possible). However, the 
algorithm is extremely powerful in practice: for instance, of the 652 635 906 combi- 
natorially distinct triangulations of the 3-sphere with 3 < n < 10 tetrahedra [10] . 
this simplification algorithm quickly reduces all but 26 of them. 

Having a powerful simplification algorithm is essential in computational 3-ma- 
nifold topology, for several reasons: 

• Many significant algorithms have running times that are exponential or 
worse in the number of tetrahedra, and so reducing this number is of great 
practical importance. 

• Good simplification may allow us to avoid these expensive algorithms en- 
tirely. For instance, in 3-sphere recognition we could simplify the triangu- 
lation and then compare it to the known trivial (< 2-tetrahedron) triangu- 
lations of S^: if they match then the expensive steps involving normal and 
almost normal surfaces can be avoided altogether |13) . 

• More generally, simplification offers good heuristics for solving the home- 
omorphism problem. Simplified triangulations are often well-structured, 
which makes Regina's combinatorial recognition routines more likely to 
identify the underlying 3-manifold (see Section |4l). Moreover, if two trian- 
gulations T and 7"' represent the same 3-manifold, then after simplification 
it often becomes relatively easy to convert one into the other using Pachner 
moves (i.e., bistellar fiips) [13j . 

Most of the individual local moves that we describe in Section 12.21 are well- 
known, and were implemented by Letscher in his early software package Normal. 
Moreover, Matveev [43j [46] and Martelli and Petronio 31] perform similar moves 
in the dual setting of special spines, and the author has described several in the 
context of census enumeration [4]. We therefore outline these moves very briefly, 
although in our general setting there are prerequisites for preserving the underlying 
3-manifold that earlier authors have not required (either because they worked in 
the more flexible dual setting of spines, or they were simply proving criteria for 
minimality and/or irreducibility) . We pay particular attention to the edge collapse 
move, whose prerequisites are complex and require non-trivial algorithms to test 
quickly. The combined simplification algorithm, whose underlying heuristics have 
been fine-tuned over many years, is described here in Section [2.31 for the first time. 

In Section [2.51 we conclude our discussion on simplification with a new technol- 
ogy that will soon make its way into Regina: simplification by breadth-first search 
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through the Pachner graph. Unlike the algorithm of Section [2. 31 this is not polyno- 
mial-time; however, it is found to be extremely effective in handling the (very rare) 
pathological triangulations that cannot be simplified otherwise. 

2.1. Generalised triangulations. A generalised 3-manifold triangulation (or just 
a triangulation from here on) is formed from n tetrahedra by affinely identifying 
(or "gluing") some or all of their 4n faces in pairs. A face is allowed to be identified 
with another face of the same tetrahedron. It is possible that several edges of a 
single tetrahedron may be identified together as a consequence of the face gluings, 
and likewise for vertices. It is common to work with one-vertex triangulations, in 
which all vertices of all tetrahedra become identified as a single point. 

Generalised triangulations can produce rich constructions with very few tetra- 
hedra. For instance, Matveev obtains 103 042 distinct closed orientable irreducible 
3- manifolds from just < 13 tetrahedra |44| . including representatives from all of 
Thurston's geometries except for 5^ x R (which cannot yield a manifold of this 
type [53]). 

Figure [T] illustrates a two-tetrahedron triangulation of the real projective space 
M.P'^ . The two tetrahedra are labelled and 1, and the four vertices of each tetra- 
hedron are labelled 0, 1, 2 and 3. Faces 012 and 013 of tetrahedron are joined 
directly to faces 012 and 013 of tetrahedron 1, creating a solid ball; then faces 023 
and 123 of tetrahedron are joined to faces 132 and 032 of tetrahedron 1, effectively 
gluing the top of the ball to the bottom of the ball with a 180° twist. 




Figure 1 . A triangulation of the real projective space RP'^ 

All of this information can be encoded in a table of face gluings, which is how 
Regina represents triangulations: 



Tetrahedron 


Face 012 


Face 013 


Face 023 


Face 123 





1 (012) 


1 (013) 


1 (132) 


1 (032) 


1 


(012) 


(013) 


(132) 


(032) 



Consider the cell in the row for tetrahedron t and the column for face ahc. If 
this cell contains u (xyz), this indicates that face abc of tetrahedron t is identified 
with face xyz of tetrahedron u, using the affine gluing that maps vertices a, b and 
c of tetrahedron t to vertices x, y and z of tetrahedron u respectively. 

For any vertex V of a triangulation, the link of V is defined as the frontier of a 
small regular neighbourhood of V. This mirrors the traditional concept of a link 
in a simplicial complex, but is modified to support the generalised triangulations 
that we use in Regina. 

A closed triangulation is one that represents a closed 3-manifold (like the example 
above): every tetrahedron face must be glued to a partner, and every vertex link 
must be a 2-sphere. A bounded triangulation is one that represents a 3-manifold 
with boundary: some tetrahedron faces are not glued to anything (together these 
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form the boundary of the manifold) , and every vertex link must be a 2-sphere or a 
disc. In either case, it is important that no edge be identified with itself in reverse 
as a consequence of the face gluings. 

Regina can also work with ideal triangulations, such as Thurston's famous two- 
tetrahedron triangulation of the figure eight knot complement [57] : these are trian- 
gulations in which vertex links can be higher-genus closed surfaces. For simplicity, 
in this paper we focus our attention on closed and bounded triangulations only, 
though most of our results apply equally well to ideal triangulations also. 

Users can enter tetrahedron gluings directly into Regina, or create triangulations 
in other ways; importing from SnapPea |59: or other file formats; building "pre- 
packaged" constructions such as layered lens spaces or Seifert fibred spaces; entering 
dehydration strings [20^ or the more flexible isomorphism signatures [13! (short 
pieces of text that completely encode a triangulation); or accessing large ready- 
made censuses that hold tens of thousands of triangulations of various types. 

2.2. Individual simplification moves. In this section we outline several indi- 
vidual local moves on triangulations. Following this, in Section 12.31 we piece these 
moves together to build Regina's full simplification algorithm. 

We describe each move in full generality, as applied to either closed or bounded 
triangulations, and without any restrictions such as orientability or irreducibility. 
For each move we give sufficient conditions under which the move is "safe", i.e., 
does not change the underlying 3-manifold. 

To keep the exposition as short as possible, we simply state these conditions 
without proof. However, the proofs are simple: the key idea is that, throughout 
the intermediate stages of each move, we never crush an edge, flatten a bigon or 
flatten a triangular pillow whose bounding vertices, edges or faces respectively are 
either identifled or both in the boundary. For detailed examples of these types of 
arguments, see (for instance) the proof of Lemma 3.7 in [Bj. 

2.2.1. Pachner-type moves. Our flrst moves are simple combinations of the well- 
known Pachner moves [491 . also known as bistellar flips. All of these moves, 
when performed locally within some larger triangulation, preserve the underlying 
3-manifold with no special preconditions required. 

Definition 2.1. Consider a triangular bipyramid: this can be triangulated with 
(i) two distinct tetrahedra joined along an internal face, or with (ii) three distinct 



tetrahedra joined an internal degree three edge, as illustrated in Figure |2 (a) A 2-3 



Pachner move replaces (i) with (ii), and a 3-2 Pachner move replaces (ii) with (i) 






(a) The 2-3 and 3-2 Pachner moves (b) The "aggregate" 4-4 move 

Figure 2. Three Pachner-type moves 
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Consider a square bipyramid (i.e., an octahedron). This can be triangulated with 
four distinct tetrahedra joined along an internal degree four edge; moreover, there 
are three ways of doing this (since the internal edge could follow any of the three 
main diagonals of the octahedron). A 4-4 move replaces one such triangulation 
with another, as illustrated in Figure [2 (b)[ 

There are two additional Pachner moves: the 1-4 move and the 4-1 move. Regina 
does not use either of these: the 1-4 move complicates the triangulation more than 
is necessary (since it introduces a new vertex, which is never needed to simplify 
a triangulation |42[ I43j). and the 4-1 move is a special case of an edge collapse 
(described later in this section). The 4-4 move is not a Pachner move, but it can 
be expressed as an "aggregate" of a 2-3 move followed by a 3-2 move. 

2.2.2. Moves around low-degree edges and vertices. It is well-known that (under the 
right conditions) minimal triangulations cannot have low-degree edges or vertices. 
One typically proves this using local moves around such edges or vertices that either 
reduce the number of tetrahedra or break some prior assumption on the manifold. 

Here we outline three such moves, which Regina uses in a more general setting 
to simplify arbitrary triangulations. Unlike the earlier Pachner- type moves, these 
moves might change the underlying 3-manifold; after describing the moves, we list 
sufficient conditions under which they are "safe" . 

Definition 2.2. A 2-0 vertex move operates on a triangular pillow formed from 
two distinct tetrahedra surrounding an internal degree two vertex, as illustrated 
in Figure 3(a)[ and flattens this pillow to a single face. A 2-0 edge move operates 



on a bigon bipyramid formed from two distinct tetrahedra surrounding an internal 
degree two edge, as illustrated in Figure [3(b)[ and flattens this to a pair of faces. 




(a) The 2-0 vertex move 






(b) The 2-0 edge move 



(c) The 2-1 edge move 



Figure 3. Local moves around low-degree edges and vertices 

Consider a tetrahedron A, two of whose faces are folded together around an 
internal degree one edge, and let A' be some distinct adjacent tetrahedron as il- 



lustrated in Figure 3(c) (for clarity, the vertices of both tetrahedra are marked in 
bold). A 2-1 edge move flattens the two uppermost faces of A' together and retri- 
angulates the remaining region with a single tetrahedron to yield a new degree one 
edge, as shown in the illustration. 

To ensure that these moves do not change the underlying 3-manifold, the follow- 
ing conditions are sufficient: 

• For the 2-0 vertex move, the two faces that bound the pillow must be 
distinct (i.e., not identified) and not both simultaneously in the boundary. 
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• For the 2-0 edge move, let g and h denote the two edges that we flatten 
together, as marked in Figure |3(b)[ These edges must be distinct and not 
both boundary. Likewise, on each side of g and h, the two faces that we 
flatten together must be distinct and not both boundary. Finally, although 
we may identify a face on one side of g and h with a face on the other, we 
do not allow all four faces to be identified in pairs, and we do not allow 
two of them to be identified if the other two are both in the boundary. 

• For the 2-1 edge move, the two edges of A' that we flatten together must 
be distinct and not both boundary, and likewise the two faces of A' that 
we flatten together must be distinct and not both boundary. Note that the 
second constraint comes "for free" as a corollary of the first. 

2.2.3. Moves on the boundary. Our next moves apply only to bounded triangula- 
tions. As before, they have the potential to change the underlying 3-manifold, and 
after presenting the moves we give sufficient conditions under which they do not. 

Definition 2.3. An book opening move operates on a face with precisely two of its 



three edges in the boundary of the triangulation, as illustrated in Figure 4(a) and 
"unfolds" the two tetrahedra on either side (which need not be distinct) to create 
two new boundary faces. A book closing move is the inverse move: it operates on 
two distinct adjacent faces in the boundary and "folds" them together so that they 
become identified as a single internal face. 



(a) The book opening and closing moves (b) One type of boundary shelling move 

Figure 4. Moves on the boundary of a 3-mamfold triangulation 

A boundary shelling move operates on a tetrahedron A that has precisely one, 
two or three faces in the boundary, and simply removes the tetrahedron from the 
triangulation (effectively "plucking it off" the boundary). This move is illustrated 
in Figure [4(b)] for the case where A has two faces in the boundary. 

The book opening move will always preserve the underlying 3-manifold. The 
book closing move will preserve the 3-manifold if (i) the two boundary faces in 
question are not the only faces in that boundary component, and (ii) the two 
vertices opposite the common edge between these faces are not already identified. 
The behaviour of the boundary shelling moves depends on the number of faces of 
the tetrahedron A that lie in the boundary: 

• If A has three boundary faces, removing A always preserves the 3-manifold. 

• If A has two boundary faces, we preserve the 3-manifold if (i) the one edge 
of A not on these faces is internal to the triangulation, and (ii) the two 
remaining (non-boundary) faces of A are not identified. 

• If A has one boundary face, we preserve the 3-manifold if (i) the one vertex 
of A not on this face is internal to the triangulation, and (ii) no two of the 
three remaining (non-boundary) edges of A are identified. 
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2.2.4. Collapsing edges. Our final move is the most powerful: it collapses an edge 
of the triangulation to a single point, and if the edge has high degree then it can 
eliminate many tetrahedra. The conditions for preserving the underlying 3-manifold 
are complex, both to describe and to test algorithmically. The details are as follows. 

Definition 2.4. An edge collapse operates on an edge e of the triangulation that 
joins two distinct vertices. It crushes edge e to a point and flattens every tetrahe- 
dron containing e to a face, as illustrated in Figure [S] 



Suppose the edge e has degree d: we denote the d tetrahedra that contain it by 
Ai, . . . , Ad, and we denote the two endpoints of e by ^ and W. 

If e is an internal edge of the triangulation (i.e., its endpoints may lie on the 
boundary but its relative interior must not), then the following conditions are suf- 
ficient to preserve the underlying 3-manifold. Section 12.41 gives details on how we 
can test these conditions efficiently. 

(1) The tetrahedra Ai, . . . , A^ must all be distinct. 

(2) The two endpoints V and W (which we already know to be distinct) must 
not both be in the boundary. 

(3) Denote the d "upper" edges in the diagram that touch V by gi,...,gd, 
and denote the corresponding "lower" edges that touch W by hi,...,hd 
(so each gi, hi pair will be merged together by the edge collapse). We form 
a multigraph F (allowing loops and/or multiple edges) as follows: 

• each distinct edge of the triangulation becomes a node of F; 

• for each i = l,...,d, we add an arc of F between the two nodes 
corresponding to edges gi and hf, 

• we add an extra node d to represent the boundary, and add an arc 
from d to every node that represents a boundary edge. 

Then this multigraph F must not contain any cycles. 

(4) In a similar way, we build a multigraph whose nodes represent faces of the 
triangulation, whose arcs join corresponding "upper" and "lower" faces that 
touch V and W respectively, and with an extra boundary node connected to 
every boundary face. Again, this multigraph must not contain any cycles. 

In essence, condition ([3]) ensures that we never flatten a chain of bigons whose 
outermost edges are both identified or both boundary, and condition ^ ensures 
that we never flatten a chain of triangular pillows whose outermost faces are both 
identified or both boundary. 



i 



V 




Figure 5. Collapsing an edge of a triangulation 
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If the edge e lies in the boundary of the triangulation then we may still be able to 
perform the move: the sufficient conditions for preserving the 3-manifold are similar 
but slightly more complex, and we refer the reader to Regina's well-documented 
source code for the details. 

2.3. The full simplification algorithm. Now that we are equipped with our 
suite of local simplification moves, we can present the full details of Regina's sim- 
plification algorithm. This algorithm is designed to be both fast and effective, in 
that order of priority, and its underlying mechanics have evolved over many years 
according to what has been found to work well in practice. 

Of course, other software packages — such as SnapPea (591 and the 3-Manifold 
Recogniser |45j — have simplification algorithms of their own. To date there has 
been no comprehensive comparison between them; indeed, such a comparison would 
be difficult given the ever-present trade-off between speed and effectiveness. 

Algorithm 2.5. Given an input triangulation T, the following procedure attempts 
to reduce the number of tetrahedra in T without changing the underlying 3-manifold: 

(1) Greedily reduce the number of tetrahedra as far as possible. We do this by 
repeatedly applying the following moves, in the following order of priority, 
until no such moves are possible: 

• edge collapses; 

• low-degree edge moves (3-2 Pachner moves, or 2-0 or 2-1 edge moves); 

• low-degree vertex moves (2-0 vertex moves); 

• boundary shelling moves. 

(2) Make up to 5R successive random 4-4 moves, where R is the maximum 
number of available 4-4 moves that could be made from any single triangu- 
lation obtained during this particular iteration of step (0). If we ever reach 
a triangulation from which we can greedily reduce the number of tetrahedra 
(as defined above) then return immediately to step (Qp. 

(3) // the triangulation has boundary, then perform book opening moves until 
no more are possible. If this enables us to collapse an edge then do so and 
return to step (QJj. Otherwise undo the book openings and continue. 

(4) // a book closing move can be performed, then do it and return to step 
Otherwise terminate the algorithm. 

The greedy reduction in step ([T]) prioritises edge collapses, because these can 
remove many tetrahedra at once, and because we typically aim for a one-vertex 
triangulation. In step ([3]), the coefficient 5 is chosen somewhat arbitrarily; note 
also that the quantity R might increase as this step progresses. The book openings 
in step aim to increase the number of vertices without adding new tetrahedra, 
in the hope that an edge collapse becomes possible. The book closures in step (|4]) 
aim to leave us with the smallest boundary possible, which becomes advantageous 
during other expensive algorithms (such as normal surface enumeration). 

Note that we never explicitly increase the number of tetrahedra (e.g., we never 
perform an explicit 2-3 Pachner move). This is only possible because we have a 
large suite of moves available: if we reformulate our algorithm in terms of Pachner 
moves alone then almost every step would require both 2-3 and 3-2 moves. 

One of the more prominent simplification techniques that we do not use is the 
O-efficiency reduction of Jaco and Rubinstein [33] . This is because the best known 
algorithms for testing for O-efficiency run in worst-case exponential time. Moreover, 
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experimental observation suggests that — ignoring well-known exceptions, such as 
reducible manifolds and MP'^ — after running Algorithm 12.51 we typically find that 
the triangulation is already O-efficient. We discuss algorithms for O-efhciency testing 
further in Section [3.31 

2.4. Time complexity and performance. Simplification is one of the most com- 
monly used "large-scale" routines in Regina's codebase, and it is imperative that 
it runs quickly. In this section we analyse the running time, which includes a dis- 
cussion of key implementation details for the edge collapse. All time complexities 
are based on the word RAM model of computation (so, for instance, adding two 
integers is considered 0(1) time as long as the integers only require logn bits). 

Throughout this discussion, we let n denote the number of tetrahedra in the 
triangulation. Some key points to note: 

• Adding a new tetrahedron takes 0(1) time. However, deleting a tetrahedron 
takes 0{n) time for Regina since we must reindex the tetrahedra that re- 
main (i.e., if we delete tetrahedron A.^ then we must rename Aj to Aj_i for 
all j>i)n Deleting many tetrahedra can be done in combined 0(n) time, 
since if we are careful in our implementation then each leftover tetrahedron 
only needs to be reindexed at most once. 

• Computing the skeleton of the triangulation (i.e., identifying and indexing 
the vertices, edges, faces and boundary components of the triangulation, 
and linking these to and from the corresponding tetrahedra) can be done 
in 0{n) time using standard depth-first search techniques. 

Theorem 2.6. Alaorithm \2.5\ (the full simplification algorithm) runs in Oiji^ logn) 
time, where n is the number of tetrahedra in the input triangulation. 

The worst culprit in raising the time complexity is the edge collapse move, and 
specifically, testing its sufficient conditions. In practice running times are much 
faster than quartic, and with more delicacy we could bring down the theoretical 
time complexity to reflect this; we return to such issues after the proof. 

Proof. First, we observe that for every move type except for the edge collapse, we 
can test the sufficient conditions in 0(1) time and perform the move in 0{n) time 
(where the dominating factor is deleting tetrahedra and rebuilding the skeleton, as 
outlined above). 

For the edge collapse move, performing the move takes 0(n) time but testing 
the sufficient conditions is a little slower. Recall that we must build a multigraph F 
and ensure that it contains no cycles. We can do this by adding one arc at a time, 
and tracking connected components: if an arc joins two distinct components then 
we merge them into a single component, and if an arc joins some component with 
itself then we obtain a cycle and the sufficient conditions fail. 

To track connected components, we use the well-known union-find data structure 
[22) . With union- find, the operations of (i) identifying which component a node 
belongs to and (ii) merging two components together each take O(logn) time. 
Therefore testing the multigraph F for cycles takes O(nlogn) time overall, and 
testing sufficient conditions for an edge collapse likewise becomes O(nlogn). 

^This is Regina's own implementation constraint: many moves could be 0(1) if we ignored 
Regina's need to consecutively index tetrahedra and other objects. Either way, however, the time 
complexity for the edge collapse — and hence the full simplification algorithm — would be the same. 
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From here the running time is simple to obtain. The algorithm works through 
"stages", where in each stage we either reduce the number of tetrahedra, or we 
reduce the number of boundary faces through a book closing move. Either way, it 
is clear there can be at most 0{n) such stages in total (since the triangulation and 
its boundary cannot disappear entirely). 

Within each stage we perform 0{n) moves in total: at most 5R < 5(# edges) < 
5 ■ 6n successive 4-4 moves, and at most (# internal faces) < 2n successive book 
opening moves. This gives us a total of 0{n^) moves throughout the life of the 
algorithm, each requiring 0{n) time to perform. 

Between each pair of moves, we might test a large number of potential moves 
whose sufficient conditions ultimately fail. The number of moves that we test on 
a given triangulation is clearly 0{n) (since there are 0{n) possible "local regions" 
in which each type of move could be performed), and so throughout the entire 
algorithm with its 0{n^) moves we run a total of 0{n^) tests. In the worst case 
(edge collapses), each test could take O(nlogn) time. 

It is clear now that the total time spent performing moves is 0{n^) and the total 
time spent testing sufficient conditions is O(n^logn), yielding a running time of 
0(n'* log n) overall. □ 

Some further remarks on this running time: 

• We noted earlier that in practice running times are faster than 0(n'* logn). 
This is because the powerful edge collapse moves typically reach the fewest 
possible number of vertices very quickly (i.e., one vertex for a closed man- 
ifold, or else one vertex for each boundary component). Once we achieve 
this, testing sufficient conditions for an edge collapse move becomes 0(1) 
time, which eliminates an n log n factor from our running time. 

• In theory, we can remove a factor of n as follows: once greedy simplification 
fails and we move on to steps ^ and ([3]) , we only test for new simplification 
moves in the immediate neighbourhood of the last 4-4 or book opening 
move. The implementation becomes more subtle and the bookkeeping more 
complex, and this is planned for future versions of Regina. 

• Finally, we note that the logn factor can be stripped down to "almost 
constant" : essentially an inverse of the Ackermann function, and < 4 in all 
conceivable situations. We can do this by applying the path compression 
optimisation to the union-find data structure; see [22], 1551 for details. 

We finish this section with a practical demonstration. Let T be the 23-tetra- 
hedron triangulation of the Weber-Seifert dodecahedral space presented in [T5] , and 
let S be an arbitrarily-chosen vertex normal surface of the highest possible genus 
(here we choose vertex surface #1733, which is an orientable genus 16 surface). If 
we cut 7" open along the surface S using Regina's cutAlongO procedure, we obtain 
a (disconnected) bounded triangulation T' with 1990 tetrahedra. The simplification 
algorithm reduces this to 135 tetrahedra in roughly 1.9 seconds (as measured on a 
2.93 GHz Intel Core i7). 

If we perform a barycentric subdivision on 7"', we obtain a new triangulation T" 
with 47 760 tetrahedra (and a large number of vertices, which means a large number 
of expensive edge collapse moves). Again the simplification algorithm reduces this 
to 135 tetrahedra, but this time takes 780 seconds, suggesting (for this arbitrary 
example) only a quadratic — not quartic — growth rate in n. 
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2.5. Exhaustive simplification via the Pachner graph. We finish our dis- 
cussion on simplification with a new technology: breadth-first search through the 
Pachner graph. The key idea is, instead of using greedy simplification heuristics, 
to try all possible sequences of Pachner moves (up to a user-specified limit). The 
result is a much slower, but also much stronger, simplification algorithm. This algo- 
rithm is based on ideas from the large-scale experimental study of Pachner graphs 
described in |13l . and the code will soon be merged into Regina's main source tree. 

Definition 2.7. For a closed 3-manifold triangulation Ai, the restricted Pachner 
graph .^i{A4) is the infinite graph whose nodes correspond to isomorphism classes 
of one-vertex triangulations of M (where by isomorphism we mean a relabelling 
of tetrahedra and/or their vertices), and where two nodes are joined by an arc if 
there is a 2-3 or 3-2 Pachner move between the corresponding triangulations. The 
nodes of l3^i{A4) are partitioned into finite levels 1, 2, . . . according to the number 
of tetrahedra in the corresponding triangulations. 

The basic idea is as follows. By a result of Matveev if we exclude level 1 

then the restricted Pachner graph ^l{^A) is always connected, and so we should 
be able to simplify a non-minimal triangulation of A4 by finding a path through 
^i(A4) from the corresponding node to some other node at a lower level. In detail: 

Algorithm 2.8. Given a one-vertex triangulation T with n tetrahedra representing 
a closed 3-manifold M, as well as a user-defined "height parameter" h, the following 
procedure attempts to find a triangulation of M with < n tetrahedra: 

(1) Conduct a breadth-first search through ^i{Ai) starting from the node rep- 
resenting T, but restrict this search to consider only nodes at levels < n + h. 

(2) // we ever reach a node at level < n, this yields a simpler triangulation 
(which we try to simplify further with the fast Algorithm \ 2.5\) . Otherwise 
we advise the user to try again with a larger h (if they can afford to do so). 

The height parameter is needed because is infinite, and because even 

the individual levels grow extremely quickly (for M = the growth rate is at least 
exponential, and it is open as to whether it is super-exponential [2]). Extremely 
small height parameters work very well in practice: for M = there arc no known 
cases for which h = 2 will not suffice |12j . 

Because of the enormous number of nodes involved, a careful implementation 
of Algorithm 12.81 is vital. We cannot afford to build even an entire single level of 
^i{M) beforehand; instead we construct nodes as reach them, and cache them 
using their isomorphism signatures (polynomial-time computable strings that also 
manage isomorphism testing [12]). The algorithm lends itself well to parallelisation 
(using multithreading with shared memory, not large-scale distributed processing) . 

Unlike the earlier Algorithm [231 this new Algorithm l2.8l is certainly not polyno- 
mial-time. Even if we are able to simplify the triangulation using k Pachner moves 
for small fc, we still test a total of 0{{n-\-h)'') potential moves; moreover, the growth 
rate of k itself is unknown, and so even the exponent could become exponential. 
See |12j for experimental measurements of k for a range of different 3-manifolds. 

In practice, with height parameter h = 2, Algorithm 12.81 easily simplifies the 
26 "pathological" triangulations of with < 10 tetrahedra that the faster Algo- 
rithm [2?5] could not. For these cases we need 5 < fc < 9 moves, with CPU times 
ranging from 0.8 to 14 seconds. This is quite slow for just n < 10 tetrahedra, and 
highlights the ever-present trade-off between speed and effectiveness. 
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3. Normal and almost normal surfaces 

One of Regina's core strengths is its ability to enumerate and work with normal 
and almost normal surfaces. A normal surface in a 3-manifold triangulation T is a 
properly embedded surface in T that meets each tetrahedron in a (possibly empty) 
collection of disjoint curvilinear triangles and/or quadrilaterals, as illustrated in 



Figure 6(a) An octagonal almost normal surface is defined in the same way, but 
also requires that exactly one tetrahedron contains exactly one additional octagonal 
piece, as illustrated in Figure |6(b)[ 



(a) Triangles and quadrilaterals (b) An octagonal piece 

Figure 6. A tetrahedron intersecting a normal or almost normal surface 

Normal surfaces are a powerful tool for high-level recognition and decomposition 
algorithms; prominent examples include unknot recognition [26], connected sum 
decomposition [33], and testing for incompressible surfaces [3T]. Almost normal 
surfaces were introduced by Rubinstein, and play a central role in algorithms such 
as 3-sphere recognition [51], Heegaard genus [37], and recognising small Seifert 
fibred spaces [52) . Rubinstein originally defined almost normal surfaces to include 
either a single octagonal piece or a single tube piece, but Thompson later showed 
that for 3-sphere recognition, only octagons need to be considered [5S] . 

We begin in Section 13.11 with a very brief overview of the necessary concepts 
from normal surface theory; for more context we refer the reader to [27) . 

In Section 13.21 we discuss the all-important problem of enumerating vertex and 
fundamental normal surfaces, and introduce a new trie-based optimisation to al- 
leviate the most severe bottlenecks in the enumeration algorithm. We follow in 
Section 13.31 with a brief discussion of 0-efficiency and the important problem of lo- 
cating normal spheres, describing the rationale behind Regina's choice of algorithm, 
and explaining why other well-known options are not effective. 

Section 13.41 outlines Regina's current implementations of 3-sphere recognition, 
3-ball recognition and connected sum decomposition. Although the key ideas are 
already known, the implementations have evolved to the point where all three algo- 
rithms are now surprisingly simple, and we present them here as a useful reference 
in a modern algorithmic form that is "ready for implementation" . 

We finish in Section [3.51 with a brief discussion of tree traversal algorithms, a new 
technology soon to appear in Regina based on backtracking and linear program- 
ming, and with enormous potential for improving performance on large problems. 

Beyond the algorithms described here, Regina offers many ways to analyse nor- 
mal surfaces, both "at a glance" and in detail. It supports the complex operation 
of cutting a triangulation open along a normal surface and retriangulating, and it 
supports the Jaco-Rubinstein operation of crushing a surface to a point [33 (which 
may introduce additional changes in topology). 
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3.1. Preliminaries from normal surface theory. In an 7i-tetrahedron trian- 
gulation 7~, normal surfaces correspond to integer vectors in a cone of the form 
{x G R^" \Ax = 0, X > 0}, where the matrix A of matching equations is de- 
rived from T. The 7n coordinates are grouped into 4n triangle coordinates, which 
count the triangles at each corner of each tetrahedron, and 3n quadrilateral coor- 
dinates, which count the quadrilaterals passing through each tetrahedron in each 
of the three possible directions. Such vectors must also satisfy the quadrilateral 
constraints, which require that at most one quadrilateral coordinate within each 
tetrahedron can be non-zero. These constraints map out a (typically non-convex) 
union of faces of the cone above. 

An important observation is that non-trivial connected normal surfaces can be 
reconstructed from their 3n quadrilateral coordinates alone ,58 . We can therefore 
identify such surfaces with integer points in a smaller-dimensional cone of the form 
{x e M'^" I Bx = 0, X > 0}. We refer to M''" and R'^" as working in standard 
coordinates and quadrilateral coordinates respectively. 

A normal surface is called a (standard or quadrilateral) vertex surface if its 
vector in (standard or quadrilateral) coordinates lies on an extreme ray of the cor- 
responding cone, and it is called a (standard or quadrilateral) fundamental surface 
if its vector lies in the Hilbert basis of the cone. The quadrilateral vertex and 
fundamental surfaces are typically a strict subset of their standard counterparts. 

Throughout this paper, we use the phrase almost normal surface to refer ex- 
clusively to the case where the extra piece is an octagon (not a tube). For al- 
most normal surfaces we introduce three additional octagon coordinates for each 
tetrahedron, yielding a cone in standard almost normal coordinates of the form 
{x S R^°" I Cx = 0, X > 0}. As before, non-trivial connected surfaces can be 
reconstructed from their 3n quadrilateral and 3n octagon coordinates [9], yielding 
a cone in quadrilateral- octagon coordinates of the form {x G R^" | -Dx = 0, x > 0}. 
We can likewise define vertex and fundamental surfaces in these coordinate systems. 

To finish, we make the well-known observation that Euler characteristic is a 
linear function in standard normal and almost normal coordinates |34) . though it 
is not linear in quadrilateral or quadrilateral-octagon coordinates. 

3.2. Enumeration. Many high-level algorithms are based on locating particular 
surfaces, which — if they exist — can be found as vertex normal surfaces, or for some 
more difficult algorithms, fundamental normal surfaces. Regina comes with heavily 
optimised algorithms for enumerating all vertex normal surfaces [8] or fundamental 
normal surfaces [H] in a triangulation, in all of the coordinate systems listed above. 

Here we focus on the vertex enumeration algorithm, which is based on the double 
description method for enumerating extreme rays of polyhedral cones [231 27] . We 
outline the double description method very briefly, and then introduce a new trie- 
based optimisation that yields significant improvements in its running time. 

In brief, the double description method enumerates the extreme rays of the cone 
{x e R'^ I Ax = 0, X > 0} by constructing a series of cones Co,Ci, . . ., where each 
Ci is defined only using the first i rows of A. The initial cone Co is simply the 
non-negative orthant, with extreme rays defined by the d unit vectors, and each 
subsequent cone Ci is obtained inductively from d^i by intersecting with a new 
hyperplane Hi. The extreme rays of Ci are obtained from (i) extreme rays of Ci-i 
that lie on Hi] and (ii) convex combinations of pairs of adjacent extreme rays of 
Ci-i that lie on either side of Hi. 
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There are significant optimisations that can be apphed to the double description 
method in the context of normal surface theory; see [71 18] for details. However, a 
major problem remains: the intermediate cones Ci can have a great many extreme 
rays (the well-known "combinatorial explosion" in the double description method) , 
and the combinatorial algorithn:0 for identifying all pairs of adjacent extreme rays 
of Ci is cubic in the total number of extreme rays. 

This cubic procedure stands out as the most severe bottleneck in the algorithm. 
The adjacency test is simple: two extreme rays xi, X2 of Ci are adjacent if and only 
if there is no other extreme ray z of Ci for which, whenever the ith coordinates 
of Xi and X2 are both zero, then the ith coordinate of z is zero also. The cubic 
algorithm essentially just iterates through all possibilities for xi, X2 and z. 

We improve this procedure by storing the extreme rays of d in a trie (also 
known as a radix tree) |54j . In our case this is a binary tree of depth d, illustrated 
in Figure [71 where at depth i the left and right branches contain all extreme rays 
for which the {i + l)th coordinate is zero and non-zero respectively^ The extreme 
rays themselves correspond to leaves of the tree at depth d. We only store those 
portions of the tree that contain extreme rays as descendants, so the total number 
of nodes is 0{d ■ # extreme rays ) and not 2"^+^ - 1. 

Legend 

Unknown 
Zero 
+ Positive 

+,+, ■ 

/ 

2 4, 1, ^ Extreme rays 

Figure 7. An example in M'' of storing extreme rays in a trie 

To identify all pairs of adjacent extreme rays of Ci, we first insert all extreme 
rays into the trie: each insertion takes 0{d) time (just follow a path down and 
create new nodes as needed). We then iterate through all pairs of extreme rays 
xi,X2, and to test adjacency we walk through the trie by (i) starting at the root, 
and (ii) whenever we reach a node at level i, if the (z -I- l)th coordinates of xi and 
X2 are both zero then we follow the left branch, and otherwise we follow both left 
and right branches in a depth-first manner. We declare xi and X2 to be adjacent 
if and only if we do not locate some other extreme ray z ^ xi , X2 during our walk. 

We can optimise this trie further: 

• At each node, we store the number of extreme rays in the corresponding 
subtree. This allows us to avoid the "false positives" xi and X2: if we are 
in a subtree containing one or both of Xi , X2 and the number of extreme 
rays is one or two respectively, then we can backtrack immediately. 



'There is an alternative algebraic algorithm; see |24l for theoretical and practical comparisons. 
'This is well-defined because each extreme ray is of the form {Av j A > 0} for some vector v. 
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• We can "compress" the trie by storing extreme rays at the nodes corre- 
sponding to their last non-zero coordinate, instead of at depth d, a useful 
optimisation given that our extreme rays may contain many zeroes. 

The appeal of this data structure is that we are able to target our search by only 
looking at "promising" candidates for z, instead of scanning through all extreme 
rays. It is difficult to pin down the theoretical complexity of our trie-based search; 
certainly it might be exponential in d (because of the branching), but it is clearly 
no worse than 0{d ■ # extreme rays), i.e., the total number of nodes. 

In practice, it serves us very well. Consider again the 23-tetrahedron triangu- 
lation of the Weber-Seifert dodecahedral space from [TH]. With the original im- 
plementation of the double description method (including all optimisations except 
for the trie-based search) , enumerating all 698 quadrilateral vertex normal surfaces 
requires 174 minutes (measured on a 2.93 GHz Intel Core 17). The new trie-based 
algorithm reduces this to 62 minutes, cutting the running time from roughly three 
hours down to just one. 

3.3. 0-efficiency. An important problem in normal surface theory is searching for 
normal spheres: this is at the heart of Jaco and Rubinstein's O-efhciency machinery 
[33] , and features in all of the high-level algorithms listed in Section 13.41 

A closed orientable 3-manifold is 0-efficient if its only normal 2-spheres are 
the trivial vertex linking spheres (which contain only triangles, and which must 
always be present). If a triangulation is not 0-efficient, then we can use Jaco and 
Rubinstein's "destructive crushing" procedure [33] : we crush the non- trivial sphere 
to a point, and then collapse away any degenerate non-tetrahedron pieces (such as 
footballs, pillows and so on) to become edges and faces, as illustrated in Figure EJ 




Figure 8. Examples of collapsing non-tetrahedron pieces 

The result is that every tetrahedron that contains a quadrilateral from the non- 
trivial sphere will disappear entirely, and the final triangulation (which might be 
disconnected) will be closed and have strictly fewer tetrahedra than the original. 
The crushing process might introduce topological changes, but these are limited to 
pulling apart connected sums, adding new 3-sphere components, and deleting S^, 
RP"^, X and/or L31 components. The crushing process is simple to implement 
(in stark contrast to the messy procedure of cutting along a normal surface), and 
with help from first homology groups any topological changes are easy to detect. 

In order to test for 0-efficiency (and to explicitly identify a non-trivial normal 
sphere if one exists), Regina uses the following result: 

Lemma 3.1. If a closed orientable triangulation T contains a non-vertex-linking 
normal sphere, then it contains one as a quadrilateral vertex normal surface. 

Proof. This result is widely known, but (to the author's best knowledge) does not 
appear in the literature, and so we outline the simple proof here. When we convert 
to vectors in standard coordinates, any connected non- vertex-linking normal surface 
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F can be expressed as a positive rational combination of one or more quadrilateral 
vertex normal surfaces, minus zero or more vertex linking spheres. Since the Euler 
characteristic x is linear in standard coordinates and x(5'^) > 0, it follows that if 
_F is a non-trivial normal sphere then some quadrilateral vertex normal surface Q 
must have xiQ) > Oj whereupon some rational multiple of Q must be a non-trivial 
normal sphere also. □ 

At present, Regina tests for O-efHciency by enumerating all quadrilateral vertex 
normal surfaces (up to multiples) and testing each. This of course is more work 
than we need to do, since we only need to locate one non-trivial sphere. We outline 
some tempting alternatives now, and explain why Regina does not use them. 

The first alternative is that, in standard coordinates, we could restrict our poly- 
hedral cone by adding the homogeneous linear constraint x ^ ^- This has been 
tried in Regina, but yields a substantially slower algorithm. Experimentation sug- 
gests that this is because (i) we are forced to work in the higher-dimensional R^" 
instead of M^"; and (ii) the constraint x ^ slices through the cone in a way that 
creates significantly more extreme rays, exacerbating the combinatorial explosion 
in the double description method. 

The second alternative is based on linear programming. Casson and Jaco et al. 
[30] have suggested (in essence) that, for each of the 3" choices of which quadrilateral 
coordinate may be non-zero in each tetrahedron, we could solve a linear program 
(in polynomial time) to maximise x over a corresponding sub-cone in M^". This is 
a promising approach, but it has a significant problem: for "good" triangulations, 
which typically are 0-efRcient, we must attempt all 3" linear programs before we 
can terminate. That is, 3" becomes a lower bound on the running time. 

Although the best known theoretical time complexity for a full vertex enumer- 
ation is slower than this [15j , in practice it is typically much faster. For example, 
when enumerating quadrilateral vertex normal surfaces for the first 1000 triangula- 
tions in the Hodgson- Weeks closed hyperbolic census [29^ (a good source of "diffi- 
cult" manifolds for normal surface enumeration), the optimised double description 
method gives a running time that grows roughly like 1.6", as shown in Figure [9] 
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Figure 9. Growth of the running time for the double description method 
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3.4. High-level algorithms. Here we present Regina's current implementations 
of the high-level 3-sphere recognition, 3-ball recognition and connected sum decom- 
position algorithms. As noted earlier, the key ideas are already known; the purpose 
of this description is to give a useful reference for these algorithms in a modern 
"ready to implement" form. 

It should be noted that all three algorithms guarantee both correctness and ter- 
mination (i.e., they are not probabilistic in nature). The 3-sphere recognition and 
connected sum decomposition algorithms include developments from many authors 
[HI inn [Ml IMl ISSl mi , and the 3-ball recognition algorithm is a trivial modifi- 
cation of 3-sphere recognition. See [43^ for a related but non-equivalent variant of 
3-sphere recognition based on special spines. 

Algorithm 3.2 (3-sphere recognition). The following algorithm tests whether a 
given triangulation T is a triangulation of the 3-sphere. 

(1) Test whether T is closed, connected and orientable. If T fails any of these 
tests, terminate and return false. 

(2) Simplify T using Algorithm ] 2. 51 

(3) Test whether T has trivial homology. If not, terminate and return false. 

(4) Create a list C of triangulations to process, initially containing just T . 
While C is non-empty: 

• Let J\f be the next triangulation in the list C. Remove J\f from C, and 
test whether Af has a quadrilateral vertex normal .sphere F. 

— If so, then perform the J aco- Rubinstein crushing procedure on F . 
For each connected component Af' of the resulting triangulation, 
simplify Af' and add it back into the list C 

— If not, and if Af has only one vertex, then search for a quadrila- 
teral-octagon vertex almost normal sphere in Af. If none exists 
then terminate and return false. 

(5) Once there are no more triangulations in C, terminate and return true. 

The key invariant in the algorithm above is that the original 3-manifold is always 
the connected sum of all manifolds in C The homology test in step (jSj) is crucial, 
since the Jaco- Rubinstein crushing procedure could silently delete S'^ x S^, RP^ 
and/or L^ i components. 

Algorithm 3.3 (3-ball recognition). The following algorithm tests whether a given 
triangulation F is a triangulation of the 3-ball. 

(1) Test whether F is connected, orientable, has precisely one boundary com- 
ponent, and this boundary component is a 2-sphere. If F fails any of these 
tests, terminate and return false. 

(2) Simplify F using Alaorithm \2.5\ 

(3) Cone the boundary ofFto a point by attaching one new tetrahedron to each 
boundary face, and simplify again. 

(4) Run 3-sphere recognition over the final triangulation, and return the result. 

For our final algorithm, we note that by "connected sum decomposition" we mean 
a decomposition into non-trivial prime summands (i.e., no unwanted S'^ terms). 

Algorithm 3.4 (Connected sum decomposition). The following algorithm com- 
putes the connected sum decomposition of the manifold described by a given trian- 
gulation F. We assume as a precondition that F is closed, connected and orientable. 
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(1) Simplify T using Alaorithm \2.5\ 

(2) Compute the first homology of T , and let r, t2 and t^ denote the rank, 
Z2 rank and Z3 rank respectively. 

(3) Create an input list C of triangulations to process, initially containing just 
T , and an output list O of prime summands, initially empty. 

While C is non-empty: 

• Let Af be the next triangulation in the list C. Remove Af from C, and 
test whether M has a quadrilateral vertex normal sphere F. 

— If so, then perform the Jaco- Rubinstein crushing procedure on F. 
For each connected component Af' of the resulting triangulation, 
simplify Af' and add it back into the list C 

— If not, then append Af to the output list O if either (i) Af has 
non-trivial homology, or (ii) Af has only one vertex and no quad- 
rilateral-octagon vertex almost normal sphere. 

(4) Compute the first homology of each triangulation in the output list O, sum 
the ranks, Z2 ranks and Z3 ranks, and append additional copies of x , 
M.P^ and L^ i to O so that these ranks sum to r, t2 and t^, respectively. 

On termination, the output list O will contain triangulations of the (non-trivial) 
prime summands of the input manifold. 

The key invariants of this algorithm are that (i) the input manifold is always the 
connected sum of all manifolds in C and O, plus zero or more x S^, M.P^ and/or 
La^i summands; and that (ii) every output manifold in O is prime and not S^. 

3.5. Tree traversal algorithms. There have been recent interesting develop- 
ments in computational normal surface theory that could allow us to move away 
from the double description method entirely. These are algorithms based on travers- 
ing a search tree |15U17| . and they combine aspects of linear programming, polytope 
theory and data structures. 

The resulting algorithms avoid the dreaded combinatorial explosion of the double 
description method; moreover, they offer incremental output and are well-suited to 
parallelisation, progress tracking and early termination. Most importantly, experi- 
mentation suggests that they are significantly faster and less memory-hungry — even 
when run in serial — for larger and more difficult problems. The code is already up 
and running, and will be included in the next release of Regina. 

Such tree traversal algorithms can be used for either a full enumeration of vertex 
normal surfaces [15j , or to locate a single non-trivial normal or almost normal sphere 
(for 0-efficiency testing and/or 3-sphere recognition) [17]. The key idea is to build 
a search tree according to which quadrilateral coordinates are non-zero in each 
tetrahedron, and to run incremental linear programs that enforce the quadrilateral 
constraints for those tetrahedra where decisions have been made, but ignore the 
quadrilateral constraints for those tetrahedra that we have not yet processed. 

For the full enumeration of vertex normal surfaces, details of the tree traversal 
algorithm can be found in [15] . To illustrate, we return again to the 23-tetrahedron 
triangulation of the Weber-Seifert dodecahedral space: whereas the trie-based dou- 
ble description method enumerates all 698 quadrilateral vertex normal surfaces in 
62 minutes, the tree traversal algorithm does this in just 32 minutes. 

For locating just a single normal or almost normal sphere, the tree traversal 
algorithm becomes extremely powerful: it can prove that this same triangulation of 
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the Weber-Seifert dodecahedral space is 0-efficient in under 10 seconds. Note that 
there is no early termination here — the tree traversal algorithm conclusively proves 
in under 10 seconds that no non-trivial normal sphere exists. 

This latter algorithm for locating normal and almost normal spheres relies on a 
number of crucial heuristics, and full details can be found in |17) . Perhaps most 
interesting is the following experimental observation: for "typical" inputs, this 
algorithm appears to require only a linear number of linear programs; that is, the 
typical behaviour appears to be polynomial-time. One should be quick to note that 
this is in experimentation only, and that the algorithm is not polynomial-time in 
the worst case. Nevertheless, this is a very exciting computational development. 

4. Combinatorial recognition 

In this brief section we outline Regina's combinatorial recognition code. Despite 
significant advances in 3-manifold algorithms, we as a community are still a long 
way from being able to implement the full homeomorphism algorithm — even simpler 
problems such as JSJ decomposition have never been implemented, and Hakenness 
testing (which plays a key role in the homeomorphism problem) has only recently 
become practical [T^ . These are the issues that we aim to address (or rather work 
around) here. 

In addition to slower but always- correct and always-conclusive algorithms such as 
3-sphere recognition and connected sum decomposition, Regina offers a secondary 
means for identifying 3-manifolds: combinatorial recognition. The central idea is 
that we "hard-code" a large number of general constructions for infinite families of 
3-manifolds (such as Seifert fibred spaces, surface bundles and graph manifolds). 
Then, given an input triangulation T, we test whether I' follows one of these 
hard-coded constructions, and if it does, "read off" the parameters to name the 
underlying 3-manifold. 

The advantages of this technique are: 

• It is extremely fast — all of Regina's hard-coded constructions of infinite 
families can be recognised in small polynomial time. 

• It allows Regina to recognise a much larger range of 3-manifolds than would 
otherwise be practically possible. 

There are, of course, clear disadvantages: 

• Such techniques require a lot of code if we wish to recognise each construc- 
tion in its full generality: Regina currently has over 25 000 lines of source 
code devoted to combinatorial recognition alone. 

• The code is only as powerful and general as the constructions that are im- 
plemented. For instance, with a handful of exceptions, Regina's recognition 
routines do not include any hyperbolic manifolds. 

• To be recognised, a triangulation must be well- structured — an arbitrary 
triangulation of even a simple manifold such as a lens space will not be 
recognised unless it follows one of the known constructions. 

Despite these drawbacks, combinatorial recognition is enormously useful in prac- 
tice. Regina's recognition code is particularly strong for non-orientable manifolds: 
of the 366 manifolds in the < 11-tetrahedron closed non-orientable census [10) . 
Regina is able to recognise all minimal triangulations for 334 of them, and is able 
to recognise at least one minimal triangulation for 357 — that is, all but nine. 
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If Regina cannot recognise the manifold from an input triangulation T (i.e., the 
combinatorial recognition is inconclusive), then often a good strategy is to modify 
T so that the triangulation becomes more "well-structured" . This might include 
(i) simplifying T using Algorithm 12. 5[ or (ii) the slower but stronger technique of 
performing a breadth-first from T through Pachner graph until we reach a trian- 
gulation that can be recognised, following the discussion in Section [2.51 

Regina is not the only software package to employ combinatorial recognition: 
there is also the 3-Manifold Recogniser by Matveev et al., which has extremely 
powerful recognition heuristics that can recognise a much wider range of 3-manifolds 
than Regina can. See [43l |44l |45] for details. 

5. Angle structures 

In addition to normal surfaces, Regina can also enumerate and analyse angle 
structures on a triangulation T. An angle structure assigns non-negative internal 
dihedral angles to each edge of each tetrahedron of T, so that (i) opposite edges of a 
tetrahedron are assigned the same angle; (ii) all angles in a tetrahedron sum to 2tt\ 
and (iii) all angles around any internal edge of T likewise sum to 27r (see Figure ITU)) . 
Such structures are often called semi-angle structures |35) , to distinguish them from 
strict angle structures in which all angles are strictly positive. Note that, by a simple 
Euler characteristic computation, an angle structure can only exist if 7" is an ideal 
triangulation with every vertex link a torus or Klein bottle. 




Figure 10. The conditions for an angle structure on a triangulation 

Angle structures were introduced by Casson and further developed by Lackenby 
and Rivin |38i I50| , and are a simpler (but weaker) combinatorial analogue of a com- 
plete hyperbolic structure. Some angle structures are of particular interest: these 
include taut angle structure^ in which every angle is precisely or tt (representing 
"flattened" tetrahedra) |28l I39j , and veering structures which are taut angle struc- 
tures with powerful combinatorial constraints [I1I2H]. AH of these objects have an 
interesting role to play in building a complete hyperbolic structure on T [25l [28l [35] . 

Conditions (i)-(iii) above map out a polytope in R'^", where n is the number 
of tetrahedra; the vertices of this polytope are called vertex angle structures^ and 
their convex combinations generate all possible angle structures on T. For many 
years now, Regina has been able to enumerate all vertex angle structures using the 
double description method, as outlined in Section [321 Moreover, it can detect taut 
angle structures and (more recently) veering structures when they are present. 

A newer development is that Regina can enumerate only taut angle structures. 
Detecting even a single taut angle structure is NP-complete [19]; nevertheless. 



^We follow the nomenclature of Hodgson et al. I28| — these are slightly more general than the 
original taut structures of Lackenby [39) . who also adds a coorientation constraint. 
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Regina can enumerate all taut angle structures for relatively large triangulations — 
in ad-hoc experiments it can do this for 70-80 tetrahedra in a matter of minutes. 
The underlying algorithm is based on the following simple observation: 

Lemma 5.1. Every taut angle structure is also a vertex angle structure. 

Proof. Describing angle structures by vectors in M^" as outlined above, suppose 
that T = Aai + (1 — A)a2 where A S (0, 1), and where r, ai, a2 are angle structures 
with T taut. Then both ai and a2 must have dihedral angles of zero wherever r 
has a dihedral angle of zero, whereupon it follows that ai — a2 — t. □ 

Algorithm 5.2. Given an n-tetrahedron triangulation T , the following algorithm 
enumerates all taut angle structures on T. 

First, we projectivise the polytope described by conditions (i)-(iii) above. This is 
a standard construction: we add a (3n+l)i/i coordinate, embed the original polytope 
V in the hyperplane x^n+i — 1, (ind build the cone from the origin through V . This 
replaces our bounded polytope V C R'^" with a polyhedral cone C C R'^"+i of the 
form C = {x G IR'^"^^ | ^x = 0, x > 0}. The vertex angle structures in the original 
polytope V now correspond to the extreme rays of the cone C. 

We run the double description method to enumerate all extreme rays ofC. Recall 
that this inductively constructs cones Cq,Ci, . . ., where Ci is obtained from Ci-i by 
intersecting with a new hyperplane Hi, and where each extreme ray of Ci is either 
(a) an extreme ray of Ci-i that lies on Hi, or (b) the convex combination of two 
adjacent extreme rays Xi,X2 of Ci-i that lie on opposite sides of Hi. 

Here we introduce a simple but crucial optimisation: in case (b) above, we only 
consider pairs xi , X2 that together do not have positive values in more than one 
coordinate position per tetrahedron. 

This optimisation is similar to Letscher's filtering method for normal surface 
enumeration [H]. It works because, if some pair Xi,X2 fails the final condition 
above, then (by virtue of the fact that all angles are non-negative and we always 
perform convex combinations) any vertex angle structure that we eventually obtain 
from a combination of Xi and X2 must have multiple non-zero coordinates in some 
tetrahedron, and so cannot be taut. 

Finally, we note that we can improve the enumeration algorithm further, both 
for enumerating taut angle structures and enumerating all vertex angle structures, 
by employing the same trie-based optimisation that we introduce in Section 13.21 

6. Experimentation 

In this penultimate section we illustrate how Regina can be used for both small- 
scale and large-scale experimentation, in the hope that readers can use this as a 
template for beginning their own experiments. 

In addition to its graphical user interface, Regina offers a powerful scripting fa- 
cility, in which most of the C-H- classes and functions in its mathematical engine are 
made available through a dedicated Python module. Python is a popular scripting 
language that is easy to write and easy to read, and the Python module in Regina 
makes it easy to quickly prototype new algorithms, run tests over large bodies 
of census data, or perform complex tasks that would be cumbersome through a 
point-and-click interface. 

Users can access Regina's Python module in two ways: 
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• by opening a Python console from within the graphical user interface, which 
allows users to directly study or modify data in the current working file; 

• by starting the command-line program regina-python, which brings up a 
standalone Python prompt. 

Users can also run their own Python scripts directly via regina-python, embed 
scripts within data files as script packets, or write their own libraries of frequently- 
used routines that will be loaded automatically each time a Regina Python session 
starts. 

The following sample Python session constructs the triangulation of M.P^ that 
was illustrated in Section [2.11 prints its first homology group, enumerates all vertex 
normal surfaces, and then locates and prints the coordinates of the vectors that 
represent vertex normal projective planes. 

babOrosemary : "■$ regina-python 
Regina 4.93 

Software for 3-manifold topology and normal surface theory 
Copyright (c) 1999-2012, The Regina development team 
>» tri = NTriangulationO 
>» to = tri .newTetrahedronO 
>» tl = tri .newTetrahedronO 

»> tO.joinToCO, tl, NPerm4(l,0,3,2)) # Glues (123) -> 1 (032) 

»> tO.joinTod, tl, NPerm4(l,0,3,2)) # Glues (023) -> 1 (132) 

»> t0.joinTo(2, tl, NPerm4(0, 1 , 2 ,3) ) # Glues (013) -> 1 (013) 

»> t0.joinTo(3, tl, NPerm4(0 , 1 , 2 , 3) ) # Glues (012) -> 1 (012) 

>» print tri . getHomologyHlO 
Z_2 

>» s = NNormalSurfaceList. enumerate (tri, NNormalSurf aceList . STANDARD, 1) 
»> print s 

5 vertex normal surfaces (Standard normal (tri-quad) ) 
>» for i in range (s .getNumberOf Surf aces ()) : 
... if s . getSurf ace (i) . getEulerCharacteristicO == 1: 

... print s . getSurf ace(i) 

0000;010||0000;010 
0000;001||0000;001 
>» 

One of Regina's most useful facilities for experimentation is its ability to create 
census data: exhaustive lists of all 3-manifold triangulations (up to combinatorial 
isomorphism) that satisfy some given set of constraints. The census algorithms are 
heavily optimised [H [HI [TO] ) and can be run in serial on a desktop or in parallel on 
a large supercomputer. 

The simplest way for users to create their own census data is through the 
command-line tricensus tool. The following example constructs all 532 closed 
orientable 3-manifold triangulations with n = 4 tetrahedra: 

babOrosemary : "$ tricensus — tetrahedra=4 — internal — orientable — finite 

— sigs output . txt 
Starting census generation. . . 

0:1 0:0 1:0 1:1 I 0:2 0:3 2:0 2:1 I 1:2 1:3 3:0 3:1 I 2:2 2:3 3:3 3:2 
0:1 0:0 1:0 1:1 I 0:2 0:3 2:0 3:0 I 1:2 2:2 2:1 3:1 I 1:3 2:3 3:3 3:2 
0:1 0:0 1:0 1:1 I 0:2 0:3 2:0 3:0 I 1:2 3:1 3:2 3:3 I 1:3 2:1 2:2 2:3 
0:1 0:0 1:0 2:0 I 0:2 1:2 1:1 3:0 I 0:3 2:2 2:1 3:1 I 1:3 2:3 3:3 3:2 
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0:1 0:0 1:0 2:0 I 0:2 1:2 1:1 3:0 I 0:3 3:1 3:2 3:3 I 1:3 2:1 2:2 2:3 

0:1 0:0 1:0 2:0 I 0:2 2:1 2:2 3:0 I 0:3 1:1 1:2 3:1 I 1:3 2:3 3:3 3:2 

0:1 0:0 1:0 2:0 I 0:2 2:1 3:0 3:1 I 0:3 1:1 3:2 3:3 I 1:2 1:3 2:2 2:3 

1:0 1:1 1:2 2:0 I 0:0 0:1 0:2 3:0 I 0:3 3:1 3:2 3:3 I 1:3 2:1 2:2 2:3 

1:0 1:1 2:0 2:1 I 0:0 0:1 3:0 3:1 I 0:2 0:3 3:2 3:3 I 1:2 1:3 2:2 2:3 

1:0 1:1 2:0 3:0 I 0:0 0:1 2:1 3:1 I 0:2 1:2 3:2 3:3 I 0:3 1:3 2:2 2:3 
Finished. 

Total triangulations : 532 

The triangulations are converted to text-based isomorphism signatures [13j and 
written to the text file output.txt, one per line. 

We can now (as an illustration) search through this census for all two-vertex 
0-efficient triangulations of the 3-sphere: 

babOrosemary : -$ regina-python 
Regina 4.93 

Software for 3-manifold topology and normal surface theory 
Copyright (c) 1999-2012, The Regina development team 
>» f = open (' output .txt ' , 'r') 
>» sig = f.readlineO 
>» while sig: 

sig = sig[0:-l] # Strip off trailing newline 
tri = NTricOigulation.fromlsoSigCsig) 
if tri . getNumberOf VerticesO == 2: 

if tri . isZeroEf f icient and tri . isThreeSphere () : 
print sig 
sig = f.readlineO 

eLAkaccddjgjqc 
>» 

Here we see that, for rt = 4 tetrahedra, there is one and only one such triangula- 
tion. To study it in more detail, we can open up Regina's graphical user interface 
and create a new triangulation from the isomorphism signature eLAkaccddjgjqc. 

These small examples illustrate how, using the census facility and Python script- 
ing combined, Regina can be an invaluable tool for experimentation, testing con- 
jectures, and searching for pathological examples. 



7. The future of Regina 

Regina continues to enjoy active development and regular releases. The develop- 
ers are currently working towards a major version 5.0 release, which will also work 
with triangulated 4^-manifolds and normal hypersurfaces (in joint work with Ryan 
Budney). Other coming developments include richer operations on triangulated 2- 
manifolds, visualisation of vertex links in 3-manifold triangulations (by Budney and 
Samuel Churchill), and much more sophisticated algebraic machinery in 2, 3 and 4 
dimensions (by Budney). Much of this code is already running and well-tested. 

Users are encouraged to contribute code and offer feedback. For information on 
new releases, interested parties are welcome to subscribe to the low-traffic mailing 
list regina-announceOlists . sourcef orge .net. 
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